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Abstract 

We use the formalism of Geometrothermodynamics to describe chemical reactions in the context 
of equilibrium thermodynamics. Any chemical reaction in a closed system is shown to be described 
by a geodesic in a 2— dimensional manifold that can be interpreted as the equilibrium space of the 
reaction. We first show this in the particular cases of a reaction with only two species corresponding 
to either two ideal gases or two van der Waals gases. We then consider the case of a reaction with 
an arbitrary number of species. The initial equilibrium state of the geodesic is determined by the 
initial conditions of the reaction. The final equilibrium state, which follows from a thermodynamic 
analysis of the reaction, is shown to correspond to a coordinate singularity of the thermodynamic 
metric which describes the equilibrium manifold. 
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I. INTRODUCTION 



The geometric nature of the thermodynamics has been considered as an important ques- 
tion since the pioneering works of Gibbs and Caratheodory However, it was the 
development of differential geometry and Einstein's theory of General Relativity which in- 



Dranches of 
SHSj set up 



creased the interest to extend and develop the geometric approach to other 
physics. Particularly, in thermodynamics, the work of Hermann and Mrugala 
the manifold called thermodynamic phase space, where the "contact geometry" approached 
by Gibbs and Caratheodory becomes meaningful. Riemannian geometry was first introduced 



in the space of equilibrium states by Rao js], in 1945, by means of a metric whose compo- 
nents in local coordinates coincide with Fisher's information matrix. In fact, the Fisher-Rao 
metric can be considered as an element of the class of so-called Hessian metrics whose local 
components coincide with the Hessian of any arbitrarily chosen thermodynamic potential, 
■lao's original work has been followed up and extended by a number of authors (see, e.g., 
3] for a review). The proposal coined by Quevedo js] as Geometrothermodynamics (GTD) 
was essential to unify both approaches and to endow the equilibrium states manifold with 
a Legendre invariant metric. 

The importance of Legendre invariance lies in the thermodynamics itself, meaning that 
once a representation is chosen to describe the system (for instance, the internal energy 
or entropy), its Legendre transform (e.g., the Gibbs free energy or the Massieu-Planck 
potential) contains the same information as the original representation. Therefore, Legendre 
invariance should be an essential ingredient of a geometric construction. 

Basically, in the GTD approach, the thermodynamic phase space is endowed with a 
Legendre invariant metric, and its maximally integral submanifold, that inherits its metric 
structure, is identified with the space of equilibrium states. In essence, a point of this space 
corresponds to an equilibrium state and, therefore, the thermodynamic processes take place 
in the equilibrium manifold. Consequently, one expects the geometric properties of the 
equilibrium manifold to be related to the macroscopic physical properties. The details of 
this relation can be summarized in three points: 

• The curvature of the equilibrium manifold reflects the thermodynamical interaction. 

• The phase transitions correspond to curvature singularities. 
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There exists a correspondence between quasi-static thermodynamic processes and cer- 
tain geodesies of the equihbrium manifold. 



The geometrothermodynamic approach has been apphed to classical systems such as the 
ideal gas js] and the van der Waals gas [o], to more exotic systems such as black holes 10 1, 
and also in the context of relativistic cosmology to describe the evolution of our Universe 



In all the cases, in which the analysis have been performed completely, the summarized 



items have been confirmed and the results have been shown to be Legendre invariant. 

In this paper, we aim to describe the geometry behind a chemical reaction. First, we will 
consider the case of a reaction with only two species, and then we will show that results can 
be generalized to include any arbitrary finite number of species. The paper is organized as 
follows. In section |Tll we present a review of the geometrothermodynamic structures and 
the classical thermodynamics approach to chemical reactions; particularly, we analyze the 
reaction A^^^ ^ B^^^ in the context of classical thermodynamics. In sections IIIII and \IV\ 
we consider the case of a reaction with two species corresponding to ideal gases and van 
der Waals gases, respectively. In both cases, we present the thermodynamic analysis of the 
reaction to find the final state of equilibrium, and analyze the same situation from the point 
of view of GTD. It is shown that each reaction can be represented by a geodesic in the 
equilibrium manifold. Then, in section |Vl we show the applicability of GTD to a general 
chemical reaction. Finally, section I VI I is devoted to the conclusions. 



II. THE THEORY 

A. Geometrothermodynamics 

The idea behind the geometrization of a thermodynamic system is simple: to build a space 
where each point corresponds to an equilibrium state. The physics behind the equilibrium 
thermodynamics allows us to say that this space is an n-dimensional manifold with the 
dimension corresponding to the number of thermodynamic degrees of freedom of the system. 
In consequence, one needs only n independent variables to coordinatize the manifold. 

In standard equilibrium thermodynamics, to a system with n degrees of freedom it is 
possible to associate n extensive variables E"-, n intensive variables J", where the index a 
runs from 1 to n, and a thermodynamic potential $, relating them. In this context, the terms 
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extensive and intensive are general concepts that refer to the independence or dependency 
of the variables associated to a given potential. For example, for a closed simple system 
with two degrees of freedom, the independent variables are T and P, if the potential chosen 
to describe the system is G, or U and V, if the fundamental potential is 5*. Recall that 
S = S{U,V). 

Consequently, from the point of view of thermodynamics, to a system with n degrees of 
freedom we associate 2n + l variables, n of them being independent. Geometrically, this idea 
corresponds to an embedding if of an n-dimensional manifold S into a {2n + l)-dimensional 
manifold T given by 

if-.s-^r, (1) 

or, in coordinates, 

^ : {E"} {$(^"), I\E^), E"} , (2) 

where b also goes from 1 to n. The manifold T is a contact manifold . This means that T 
is endowed with a family of tangent hyperplanes (contact structure) defined by the so-called 
fundamental 1-form O that satisfies the non-integrability condition 

e A (de)" y^O . (3) 

The manifold T is called the thermoynamic phase space. It can be coordinatized by the 

2n + l variables {Z^} = {$, E"", where A = 0, 2n. Its existence is relevant because it 

results adequate to perform Legendre transformations (as in thermodynamics) as a change 

of coordinates. Formally, a Legendre transformation is a contact transformation, i.e., a 

transformation which leaves the contact structure unchanged; in coordinates 1 2'"^), it is 
□ 

defined as 

{Z^} {Z^} = {$, r} (4) 

$ = -I - SkiE'^P , E' = -P , E^ = & , P = E\ P = P , (5) 

where iU j is any disjoint decomposition of the set of indices {1, . . . , n}, and k,l = 1, . . . ,i. 

According to the Darboux theorem , the 1-form of equation ([3]) can be given in the 
coordinates {Z^} as: 

e = d$ - PdE"" (6) 



where we use Einstein's summation convention for repeated indices. It can easily be seen 
that after a Legendre transformation, the new 1-form G in coordinates {Z"^} reads 



e = d<l - ladE" . (7) 

This proves that the contact structure remains unchanged. 

On other hand, the manifold S is the maximally integral submanifold of T, and is defined 
in such a way that the properties of the thermodynamic systems are encoded in it. So far, 
we have seen that S is specified through the embedding ([1]), which is equivalent to specifying 
the fundamental equation of the system ^{E"-). The next step is to introduce the relations 
of standard equilibrium thermodynamics into the definition of the manifold. To this end, 
we demand that the embedding ([1]) satisfies the condition 

^*{Q) = , (8) 

where cp* is the pullback of (p. In coordinates, it takes the form 

^*(e) = ^* (d<|. - lAE^ = (J^ - l}j dE^ = 0. (9) 



It follows immediately that 



9$ 

^ = ^E^) and g^ = Ia. (10) 



Equations (jH]) and fllOp constitute the standard Gibbs relations of equilibrium thermody- 
namics in S, namely, 

d$ = ladE" . (11) 

The manifold S defined in this way is called the space of equilibrium states. 

In addition to the geometric description of thermodynamics in terms of a contact struc- 
ture, the GTD program promotes the contact manifold (T, 0) into a Riemannian contact 
manifold {T,Q,G), where G is a metric sharing the symmetries of 0. The most general 
metric invariant under total and partial Legendre transformations that has been found so 



far is given by 
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n 

G = e ® e + A [Z^] [(^"4)'"^' d^" ® dia] , (12) 



a=l 



where A{Z^) is an arbitrary Legendre invariant function of the coordinates and k is an 
integer. The corresponding induced metric in the space of equihbrium states is given by 

5$ \ """^^ 



9^ = ^*iG) = AiE'^)J2 



(13) 

a,b=l L ^ 

As suggested in [13(, this metric could be useful to analyze multicomponent systems, par- 
ticularly systems where chemical reactions take place. We will show in the next sections 
that, in fact, chemical reactions can be represented as geodesies of the equilibrium manifold 
described by the metric f|T3l) . 

Notice that to compute the explicit components of this metric, it is necessary to specify 
only the fundamental equation $ = $(£'"). Thus, all the geometric properties of the equi- 
librium space are determined by the fundamental equation. This is similar to the situation 
in classical thermodynamics where the fundamental equation is used to determine all the 
equations of state and thermodynamic properties of the system. 

Notice that the metric G contains the arbitrary parameter k which, however, can be 
absorbed by renaming the coordinates as dX"- = [E°-)'^'''^^dE°- and dYa = {laY^^^dla- Then, 

n 

G = e ® + A(z^) 5^ (dx'^ ® di;) . (14) 

a=l 

Furthermore, the arbitrary function A(Z"^) can be fixed by demanding invariance with re- 
spect to changes of representation, an issue which is outside of the scope of the present work 



14| . It is therefore possible to perform the entire analysis with an arbitrary function A(Z^) 
in coordinates = ($,X°,Ya); however, the physical interpretation of these coordinates 
becomes cumbersome and makes it difficult the physical interpretation of the results. For 
the sake of simplicity, we will use in this work the particular choice k = —1 and A = — 1, 
which has been shown to be useful also to describe geometrically systems like the ideal gas 
or van der Waals gas [l^. Then, in the particular case n = 2, the metric of the equilibrium 
manifold reduces to 



E^—r + E^ 



dE^ / \ dE 



dE'® dE^ . (15) 



dE^dE 



To analyze the geometric properties of the equilibrium manifold, we will consider the 
connection and the curvature. In particular, the connection is used to represent the geodesic 
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equations in a given coordinate system 

rf^E" ^„ dE^dE" 1 ,.fdgdb dgdc dgtA 

^^'-iFur^'' = '''' 

where F'^^^ are the Christoffel symbols. The solutions of these equations are the geodesies 
E'^ij), where r is an affine parameter along the trajectory. One of the main goals of the 
present work is to show that a given chemical reaction can be represented geometrically as 
a family of geodesies of the equilibrium manifold, which is determined by the fundamental 
equation of the chemical system. 
The curvature tensor is defined as 

pa bd be I pa pe pa pe ('\'7\ 

In GTD, the curvature tensor of the equilibrium manifold is expected to be a measure of 
the interaction between the components of the thermodynamic system. Furthermore, from 
the curvature tensor one can define the Ricci tensor Rab = g'^'^Racbd and the curvature scalar 
R = Notice that in the case of a two-dimensional space, the curvature tensor has only 

one independent component, say -R1212 and, therefore, the Ricci tensor and the curvature 
scalar are proportional to -Ri2i2- 

B. Thermodynamics of chemical reactions 

Consider the general chemical reaction 

a^Ai + a^A^ + . . . ^ + baBs + . . . , (18) 

in which the ai, 02,... are the stoichiometric numbers of the reactants Ai, A2,..., and 
bi, b2,... are the stoichiometric coefficients of the products Bi, B2,..., respectively. Notice 
that the species need not all to occur in the same phase. The main condition for the 
chemical-reaction equilibrium in a closed system is that 

i 

where the coefficients Ui = (— ai, — a2, 61, 62, •••) refer to the stoichiometric numbers and 
/^i = ifJ'Ai, fJ'Bi, ■■■) represents the chemical potential of the i— species. Notice that the 
condition (fT9l) holds no matter how the closed system reaches its final equilibrium state. 



During a reaction, the change of numbers of moles of the species i, Arii = rii — Uifi, where 
Ui^o is the number of moles of species i at the beginning of the reaction, is proportional 
to the stoichiometric number i/j with the extent of reaction ^ as the proportionality factor, 
Arii = Ui^. For an infinitesimal extent of reaction dC,, it holds that drii = Vid^. The fact 
that the extent of reaction can be treated as an infinitesimal quantity is essential for the 
geometric description we will present below. 

The main premise to apply classical thermodynamics in closed systems, where chemical 
reactions can occur, is that we can use thermodynamic variables to describe the system even 
if it is not in material equilibrium. It means that variables such as [/, S*, T, V, etc., are 
completely defined at any extent of reaction. 

Though f|T9l) is useful for practical situations, it does not contain information about the 
behavior of the different thermodynamic potentials from initial conditions to equilibrium. 
To obtain this information, we will use the fundamental equation of the chemical system. 
Let (^Aji^E^^) represents the fundamental equation of the species Aj. Then, the fundamental 
equation for the general reaction (fTSjl can be constructed as follows 

m\,^El^, ...) = ^aAEX) + ■■■ + ^bAE^bJ + ■■■ + ^A.MEX.EX) + ... 

+ ^A,MEX,E^^J + ... + ^a,,a,,bAEm^E'X,,E%^) + ... (20) 

Notice that in this expression we are taking into account all possible interactions between 
all the species. The only assumption is that the interaction between the species Aj and 
Ak depends on the variables and E'\^, only. It seems that this condition is not very 
restrictive in realistic situations. 

In the last subsection, we emphasized the role of the Legendre transformations from the 
geometrical point of view. Now, we will mention their importance from the thermodynamic 
point of view. 

1. It does not matter which potential is chosen to describe a particular system, 
all of them will contain the same thermodynamic information. 

2. The prediction of the final equilibrium state is made in accordance with the "extremum 
principles" ; these principles could be different for different potentials. The importance 
of the Legendre transformations is that they always can be used to find a potential in 
which the extremum principles hold for the experimental working conditions. 
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For later use, we summarize here the "extremum principles" 
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• Entropy Maximum Principle. The equilibrium value of any internal unconstrained 
parameter is such as to maximize the entropy for the given value of the total internal 
energy. 

• Energy Minimum Principle. The equilibrium value of any internal unconstrained pa- 
rameter is such as to minimize the energy for the given value of the total entropy. 

• General Minimum Principle for the Legendre Transforms of the Energy. The equilib- 
rium value of any unconstrained internal parameter is such as to minimize the Legendre 
transform of the internal energy for a constant value of the transformed variable (s). 

• General Maximum Principle for the Legendre Transforms of the Entropy. The equi- 
librium value of any unconstrained internal parameter is such as to maximize the 
Legendre transform of the entropy for a constant value of the transformed variable(s). 

Consequently, we can talk about Legendre invariance in two senses. Firstly, in the sense 
that the thermodynamic information remains conserved under a Legendre transformation 
and, secondly, in the sense that the same equilibrium value for one or several internal 
unconstrained parameters will be obtained, independently of the thermodynamic potential 

under the condition that the experimental restrictions are according to the particular 
restrictions contained in the "extremum principle" for $. 

In the following sections, we will study the chemical reaction A^^-j ^ B^^-j considering A 
and B either as ideal monoatomic gases or as van der Waals monoatomic gases. 



III. IDEAL GASES 



A. Thermodynamics 



If we consider the species as ideal gases, t 
species in the entropy representation reads 

S{U, V, n) = nsQ + uR In 



le corresponding fundamental equation for each 
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U 



V\ ( n 



(21) 
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where sq, Uq, Vq and no refer to the standard values of reference, c is a dimensionless 
constant related to the heat capacity of the ideal gas, i.e., Cy^n = cR, and R is the universal 
gas constant. 

Let us consider the particular case in which only two species take part in the reaction. 
It turns out that it is convenient to study the evolution of the total entropy in terms of the 
extent of reaction. According to Eq.( l20|) . the fundamental equation takes the form 



S {Ua, Va, Ub, Vb, = S{Ua, Va, + S{Ub, Vb, 

Ua V ( V, 



= {^Afi - i)sQ^A + {riAfi - O^ln 
+ {nB,o + i)sQ,B + {uBfl + 



0,A 



u 



B 



0,B 



CB 



Vo,A 

Vb 



Vr 



0,B 



no,A 

nB,o + ^ 
no,B 



^(CA + I) 



-(CB+I)' 



(22) 



where we neglected the interaction term for simplicity, and introduced explicitly the variable 
extent of reaction ^. Note that rij^o refers to the initial conditions and no,i refers to the 
values of the state of reference. For s.mpl.city the vah,es for n,. and V„,. are set equaUo 
one. Moreover, the values for Uo^i and so,i -which depend on the nature of each gas [17J- 
are chosen in such a way that they basically take into account the differences between the 
species. These and other values are shown in Table HI Finally, we take the temperature as 



Gas 


A 


B 


Initial moles (mol) 


1 





Heat capacity c 


3 
2 


3 
2 


Molar standard entropy (J/mol-K) 


1 


2 


Standard internal energy (J) 


1 


2 



TABLE I: Conditions for the gases A and B 



Ta = Tb = T = 300K and the volume as Va = Vb = V = 20L, experimental conditions 
that can easily be achieved. With these conditions, the fundamental equation reduces to 

>g(0 = l + e + i^(l-01n( j+R^\ni J , (23) 

where the value of the constants have been rounded to simplify the presentation. A plot of 
this function is displayed in figure [Tal 
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Since the temperature is constant and ca = cb, the total internal energy is constant. 
Thus, according to the "entropy maximum principle" , we have that the maximum of S, as 
a function of the extent of reaction, corresponds to the equilibrium condition. In this case, 
S reaches it maximum value SfatC,=^f^ 0.285. 

To illustrate the significance of Legendre invariance, we can analyze other potentials that 
are obtained from S by means of Legendre transformations. Consider, for instance, the 
Massieu potential (Helmholtz free energy) 

1 



S 



p_. 
T 



i 

Considering the assumption fl20|) with no interacting term, we obtain 



(24) 



V, Ha, ub) = 0a(/3, V, ua) + V, ub) 

where (5 = —. To calculate explicitly this function we use 

dS- 

Hi{Ui,V,ni) = -T— ^ 

Ofli 



(25) 



Pi{Ui,V,ni) 
CiRrii 



T 



dV 



/3 



In this way, the fundamental equation of the system is: 
(t){l3,V,nA,nB) = ^ ni <( so,i + -Rln 



i=A,B 



rij 



-1 



V_ 



-cR 



(26) 
(27) 
(28) 

(29) 



Taking into account the considerations discussed for the entropy representation, and the 
value P = -;^K~^, finally the fundamental equation in this representation is reduced to 



300 

m 



-11.47 + ^ + (1 -Oln 



4.58 X 10^ 



+ In 



1.62 X 10^ 



(30) 



1-e ; V e 

where again we have rounded the values of the constants for simplicity. The plot of this 
function is displayed in figure [ibl In this case, the transformed variable of the Legendre 
Transform is /3, and since it is a constant, the maximum of corresponds to the equi- 
librium state, which is reached at ^ = ^j- ^ 0.285. Thus, in agreement with the Legendre 
invariance, we obtain the same result as in the entropy representation. 
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(a) Evolution of the Entropy (b) Evolution of the Massieu 
potential potential 

FIG. 1: Behaviour of the thermodynamic potentials for the reaction A^^^ ^ Bj-^^ at T=300 

K considering A and B as ideal gases. 

The main result of this thermodynamic analysis is that a system with initial parameters 
as given in Tabled which correspond to an entropy Sin, undergoes a chemical reaction whose 
final state is the equilibrium state characterized by the reaction extent ^ ~ 0.285 and the 
entropy Sj. Notice that for < ( 6n > ^f) all the states with ^ > ^/ < 0) 
not permitted. In fact, once the system reaches the final equilibrium state at ^ = the 
chemical reaction ends, and the states characterized by a decrease of entropy are unphysical 
according to the second law of the thermodynamics. 

B. Geometrothermodynamics 

Recall that to construct the metric of the equilibrium manifold we only need the funda- 
mental equation. Under the restrictions corresponding to the chemical reaction of two ideal 
gases as described in the last subsection, the original fundamental equation (1221) reduces to 
a function that depends on two variables only, namely, 

^(^,0 = l + e + (l-0^1n(^^)+ei?ln(^^^^^j , (31) 
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dr2 ' ' U \dT 



where we used the variable U = Ua + Ub = cijiA + nB)RT to rewrite the variables Ua and 

f/fi as 

Ua= U , Ub= + ^ U. (32) 

nA,Q + ?^B,o J^yi.o + nsfi 

In the entropy representation $ = S*, we choose the independent variables as E"- = {U,^}. 
Then, the metric of the equilibrium space (fT5!) for the chemical reaction of two ideal gases 
reduces to 

ig _ dU"^ Rd^"^ 

- u2 ^2(i_^) [i_^in(2v^0 + ^ln(l-0] ' 
For the metric ( l33l) the only non- vanishing Christoffel symbols are and r|^. Then, the 
geodesic equations read 

:^^W-^f:^f^wV = o, (34) 

+ — ^ -r 77?7Y^ ( :r^(^)) =o- (^5) 

The main point now is to see whether the geodesic equations can reproduce the thermo- 
dynamic process that occur during a chemical reaction. The idea is that we use the initial 
values of the thermodynamic variables, corresponding to the initial equilibrium state of the 
chemical reaction, to identify a particular point in the equilibrium manifold. This point is 
then used as initial value to integrate the geodesic equations. The question is whether the 
solution of the geodesic equations passes through the final equilibrium state of the chemical 
reaction. 

For the particular case of ideal gases we are investigating here, we found in the last 
subsection that the thermodynamic analysis establishes the value of the extent of reaction 

^ 0.285 for the final state. The values of the initial state have been incorporated in the 
fundamental equation (|3T|) and, consequently, in the thermodynamic metric fl33|) and in the 



^2 i?-(2-30(l-i?ln^) 



geodesic equations. We now consider the "experimental" condition T = 300Ji . Then, we 
3 

get U{t) = -8.314 300 J = const., so that equation flMl) is satisfied identically. Then, we 
proceed to solve numerically the remaining geodesic equation f l35p for ^. The results are 
displayed in figures [2al and [2bl We choose as initial conditions values very close to ^(0) = 
and .^(0) = 1 (since the reaction can go in both directions), and different arbitrary initial 
"velocities" d^{Q)/dT = ^(0). The important result is that all the geodesies reach the point 
C,f ^ 0.285, independently of the initial values of ^ and ^, and all of them follow a pattern 
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in the physical region between the initial value and where the entropy increases. In fact, 
the numerical integrator finds always a "singularity" at the point C, ~ 0.285. We will see 
below that in fact this point corresponds to a coordinate singularity of the thermodynamic 
metric. 

We now test the Legendre invariance of our analysis by using the Massieu potential fl29|) . 
We insert the conditions for the present chemical reaction and obtain the corresponding 
fundamental equation. In this case, E°- = and the metric flT^ leads to 



,g_d(3' Rde 



(36) 



The only non-vanishing Christoffel symbols are F^^ and r|^. Then, the geodesic equations 
are 



re(r) ^ ^^ir)) =0. (38) 



We now fix the "experimental" condition /3(r) = K ^, with that the equation f p7|) 

300 



is satisfied identically and we proceed to solve numerically the remaining equation fl38|) for 
^. The initial conditions are the same as in the entropy representation. The results are 
displayed in figures [3a] and [3bl We see again that all the geodesies pass through the point 
^ 0.285, independently of the initial values. Moreover, in the physical region, before the 
final equilibrium state is reached, they all coincide with the geodesies shown in Figs. [2a] and 
l2b] for the analysis in the entropy representation. This shows explicitly that the analysis 
does not depend on the choice of thermodynamic potential. 

A straightforward computation shows that the curvature tensor vanishes in both repre- 
sentations and, consequently, the corresponding equilibrium space is fiat. This indicates that 
no thermodynamic interaction exists in a chemical reaction in which only non-interacting 
ideal gases are involved. 
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(a) Solution of geodesic equation ([35 
for ^(0) = 0.01 and different initial 
"velocities" ^(0) 
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(b) Solution of geodesic equation ([35 
for ^(0) = 0.99 and different initial 
"velocities" ^(0) 



FIG. 2: Behaviour of the geodesic solution in the Entropy representation, for the reaction 



A 



(g) 



Bj-g^ at T=300 K considering A and B ideal gases 



IV. VAN DER WAALS GASES 



A. Thermodynamics 



A more realistic gas is described by the van der Waals fundamental equation [16 1 



S{U, V, n) = nso + nR In 



CRT. vAn 



(39) 



where a and b are constants. For simplicity, we will consider both gases with the same 
a value, so that the coupling terms in the fundamental equation (120!) can be considered 
as vanishing. That is, this case is a simple mixture in which the interactions A — B are 
identical to the interactions A — A and B — B. Consequently, the fundamental equation of 
the mixture reads 



S{Ua, Ub, V,nA,nB) = Sa{Ua, V,nA) + Sb{Ub, V^ub) ■ 



(40) 



In the previous section, the constant temperature condition of the ideal gas reaction was 
equivalent to implying Utotai =const. and, therefore, it was straightforward to compare the 
S and representations. In fact, the conditions, under which the extremum principle is 



15 




0.8 



0.6 



4(»io/) 



0.4 



0.2 




8 10 



4 6 



8 10 



(a) Solution of geodesic equation ([37 
for ^(0) = 0.01 and different initial 
"velocities" ^(0) 



(b) Solution of geodesic equation (f3^ 
for ^(0) = 0.99 and different initial 
"velocities" ^(0) 



FIG. 3: Behaviour of the geodesic solution in the Massieu Potential representation, for the 
reaction Aj.^^ ^ B^^-, at T=300 K considering A and B ideal gases 



'(g) 



valid, were fulfilled in both representations. On the other hand, in the case of van der Waals 
gases, the constant temperature condition does not imply that the total internal energy 
is constant. Consequently, it is necessary to be cautious when applying the extremum 
principle corresponding to a constant temperature condition. To this end, we perform a 
Legendre transformation in the S"— representation to obtain the Massieu potential so that 
the temperature is an independent variable. Then, for the chemical reaction of two van der 
Waals gases we obtain the thermodynamic potential 



V, ua, ub) = (t>A{(3, V, ua) + 4>b{(3, V, ub) 

= -(3 [nA{(3,V,nA)nA - PaV + HB{(3,V,nB)nB - PbV] 



(41) 



For the sake of simplicity, we consider the same initial values as given in Table [T] for ideal 
gases, with the additional values a = 506 J Lmol'"^ and b = 0.050 Lmo/~^, which are quite 
reasonable for real gases 181]. The behavior of the potential is displayed in Fig. |H Notice 
that in this case the maximum is reached at = ~ 0.284 which is only slightly different 
from the ideal case, because the initial equilibrium state is the same in both cases and, at 
T = 300 K, the differences between ideal and van der Waals gases are very small. In the 
next section, we will find out if there exists a geodesic in the equilibrium space that connects 
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FIG. 4: Behavior of the thermodynamic potentials for the reaction A^^-j ^ B^^^ at 
T = 300K considering A and B as van der Waals gases. 



the initial equilibrium state with the final one. 



B. Geometrothermodynamics 



Considering the initial conditions given in Table HI the fundamental equation in this case 
can be obtained directly from Eq.( l4T]) as 



=l-ci? 1 + ln 



^ 1 -ci?ln2 + i?ln 



cR 



20-6^ 



(42) 



From here one can compute all the components of the metric tensor. The explicit expressions 
are quite cumbersome and cannot be written in a compact form. The curvature scalar is in 
general different from zero, indicating the presence of thermodynamic interaction. This is 
in accordance with the statistical approach to the van der Waals gas in which, as a result of 
the interaction between the particles of the gas, the corresponding Hamiltonian possesses a 
non-trivial potential term. 

All the Christoffel symbols are different from zero so that the geodesic equations contain 
all the independent terms 



2T 



d 



d 



dr 



d 



^ndr 



2n 



0, (43) 



. (44) 
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t t 



(a) Solution of geodesic equation (f44]) (b) Solution of geodesic equation (f44ll 
for ^(0) = 0.01 and different initial for ^(0) = 0.01 and different initial 

"velocities" ^(0) "velocities" ^(0) 

FIG. 5: Behaviour of the geodesic solution in the Massieu Potential representation, for the 
reaction A^^-j ^ B^^-j at T=300 K considering A and B Van der Waals gases 

This system is solved numerically with initial values for C, very close to .^(0) = and .^(0) = 1 

and arbitrary initial "velocities" ^(0). Moreover, for /3 we introduce the "experimental" 

condition in the form /3(0) = and /3(0) = 0. The results are displayed in Figs. I5al and 

300 

l5bl We see also in this case that all the geodesies reach their final equilibrium state for 
^ = ^ 0.284 which is the same value we obtained in the thermodynamic analysis of the 
last subsection. Moreover, the geodesies follow a pattern in the physical region, in the same 
way as in the case of ideal gases. We conclude that the geodesies can be interpreted also in 
this the geometric path that the chemical reaction follows until it reaches its final 

equilibrium state. 



V. GEOMETROTHERMODYNAMICS OF CHEMICAL REACTIONS 

In the previous sections, we analyzed two particular chemical reactions with only two 
species. In both cases, we found that the corresponding equilibrium manifold S reduces 
to a 2— dimensional manifold, once the conditions of the reaction are taken into account. 
The question arises whether the dimension of S increases as the number of reactants and 
products of the reaction increases. We will show in this section that GTD can handle in a 
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simple manner the chemical reaction of any arbitrary (finite) number of species in a closed 
system. 

Consider the chemical reaction of r species described by the variables Si, Ui, Vi, and 
rii {i = l,...,r). In the entropy representation, for instance, the thermodynamic proper- 
ties of each substance is determined by the fundamental equation Si = Si{Ui,Vi,ni). Each 
fundamental equation Si generates a 3— dimensional equilibrium space Si for the species. 
According to Eq. fl20l) . the fundamental equation of the chemical system $ = $(5*1, Sr) = 
^{Ui,Vi,ni) will depend on 3r variables. Then, the total equilibrium space S has 3r di- 
mensions. In general, the level of computational difficulty in geometry increases with the 
number of dimensions, so that for large r the calculations could easily be outside the reach 
of computational capability. However, we can use the conditions of the reaction to reduce 
the number of dimensions. A common condition for reactions involving gases is that the 
reaction occurs at constant volume so that the functional dependence of the fundamental 
equation can be reduce to 2r, i.e., $ = $(?7j,nj). Furthermore, using the definition of the 
extent of reaction parameter, Aui = Ui^, we can replace all the n^'s variables by ^, according 
to rii = Uifl + Vi^. Since rij^o and Vi are constants, the functional dependence of the funda- 
mental equation reduces to $ = ^(f/j,.^). Using the equations of state for each species, we 
can express each Ui in terms of t/, rii and other constants. For instance, in the case of ideal 
gases we have that 



so that the fundamental equation becomes $ = $(?7, ^). Consequently, the corresponding 
equilibrium manifold is 2- dimensional, independently of the number of species. 

In the case of more complicated fundamental equations, it is always possible to express 
each Ui in terms ofU and the extent of reaction in such a way that the resulting equilibrium 
manifold has only two dimensions. For instance, in the case of van der Waals gases, we 
obtain 




(45) 



Then, we can express each Ui as 




(46) 



U^{U,i) 



Ciirii^o + ViS) 



( 



j 



V 



a 



(47) 
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Let us now investigate the geodesic equations. To construct the metric g of the equihbrium 
manifold we can use the thermodynamic potential ^) or any other potential that can 
be obtained from $(?7, ^) by means of a Legendre transformation (basically, $(/3,^)). The 
results do not depend on the choice of $, because the geometric properties of £ in GTD are 
Legendre invariant. If we take the potential $([/, the geodesic equations can be written 
in general as 

f (-) + r"vv iU(r)\\ 2T" (±U(r)](4-((r))+T"(i((r)X^0, (48) 



^Ur) . Ti (±Ur))\ 2Ti, {^r)) (±Ur)) + 4„ (^^r))' ^ ^ (49) 



According to the Picard-Lindelof theorem 19|, given an initial value, i.e., ^(0) , ^(0), 
f/(0) and f/(0), if £ is smooth -as in the preceding cases-, the solution to this equation exists 
and is unique. In the case of a chemical reaction, the value of .^(0) is fixed by the initial 
equilibrium state of the reaction, but the value of ^(0) remains free. In the case of ideal 
gases, in which U is proportional to T, the initial values of the second variable U{0) and 
U (0) are fixed by the conditions of the chemical reaction. In the case of more general species, 
a Legendre transformation must be performed such that T becomes the second independent 
variable and the initial values are T(0) and T(0). 

We have seen in the examples above that the final equilibrium state of the reaction does 
not depend on ,^(0) and that the numerical integrator detects a "singularity" at that point. 
It turns out that this corresponds to a coordinate singularity of the metric. To see this, we 
calculate the component (^^) of the general metric ( fT5l) . and obtain 

This expression is valid for $ = S* and all the thermodynamic potentials that can be obtained 
from S by means of any Legendre transformations -except those which change the role of 
It is easily seen that as soon as the chemical-reaction equilibrium condition (1191) is satisfied, 
the denominator goes to zero and, thus, the metric is not well defined. This means that at 
the final equilibrium point the coordinates are not appropriate to describe the equilibrium 
manifold. Notice that this result is completely general since it does not depend on the 
particular reaction that determines the potential $. 
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The fact that the final equihbrium state is characterized by a coordinate singularity allows 
us to perform an analytical investigation of that particular point. Consider, for instance, 
the metric (I33p for two ideal gases in the S"— representation. The component g^^ presents a 
physical divergence when 1 — i?,ln(2v^0 + -Rln(l — = 0. The solution of this equation 

= 7= (51) 

with R = 8.314 gives C,f ~ 0.285 which is exactly the value obtained in the numerical 
investigating of the geodesic equations for this metric. 

Finally, let us mention an additional invariance property of the GTD approach. To 
reduce the number of independent variables of the fundamental equation in the case of van 
der Waals gases, we used the relationship P7|) for Ui{U,^) with V = const. However, due 
to the invariance of the metric g under changes of coordinates, we can also use the same 
relationship as Ui = UiiV,^). Then, the resulting fundamental equation becomes ^(V,^), 
once P is fixed in accordance with the conditions of the reaction. The corresponding metric 
can be computed and the geodesic equations can be integrated numerically for V and ^ with 
the same initial conditions for C,, and V{0) = 20 and V^(0) = 0. The resulting geodesies are 
exactly the same as the ones we obtained with $([/, ,^). 



VI. CONCLUSIONS 



In this work, we used the formalism of GTD to present a geometric representation of 
chemical reactions in closed systems. In GTD, all the information about a thermodynamic 
system is encoded in its equilibrium manifold determined by a metric which is invariant under 
Legendre transformations, i.e, its properties do not depend on the choice of thermodynamic 
potential. 

First, we consider the case of a chemical reaction with only two species corresponding 
either to ideal gases or to van der Waals gases. In the case of ideal gases, we found that 
the equilibrium manifold is fiat, independently of the thermodynamic potential. In GTD, a 
vanishing curvature means that there is no thermodynamic interaction. This agrees with the 
interpretation from the point of view of statistical mechanics and thermodynamics: With the 
statistical point of view, because the molecules of each ideal gas behave as "free particles" 
since the Hamiltonian contains exclusively the kinetic part; and with the thermodynamic 
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point of view, because there are no coupling terms in the fundamental equation ( !22|) . In 
the case of van der Waals gases, the curvature is different from zero, indicating the presence 
of thermodynamic interaction. This is also in accordance with the statistical approach to 
the van der Waals system, because the Hamiltonian contains a potential term which is 
responsible for the interaction between the molecules of the system. 

The thermodynamic analysis of the chemical reactions in both cases shows that the fi- 
nal equilibrium state is reached at a particular value of the extent of reaction parameter. 
The numerical analysis of the geodesies in the corresponding equilibrium manifolds provides 
exactly the same value of the extent of reaction for the final equilibrium state. It is evi- 
dent that all the geodesies follow a path pattern that ends at the same equilibrium state, 
independently of the initial values of ^(0) and ^(0). Moreover, the final equilibrium state is 
always denoted as a singularity by the numerical integrator. 

In order to understand our results in a more general fashion, we analyzed a general 
reaction in the context of GTD. We showed that using the conditions of thermodynamic 
equilibrium and the laboratory conditions of the reaction, it is always possible to reduce to 
two the number of dimensions of the equilibrium manifold. This is an interesting result that 
allows us to describe any chemical reaction as a geodesic curve on a 2— dimensional space. 
Finally, it was shown that the metric of the equilibrium manifold possesses a coordinate 
singularity exactly at that point where the condition for the reaction equilibrium is satisfied. 

The examples of chemical reactions presented in this work involve only gases. Never- 
theless, the generalization to include reactions involving solids or liquids is straightforward. 
Indeed, once the fundamental equations of the species are given, we can construct the cor- 
responding Massieu-Planck potential of the reaction which allows us to consider P and T 
as constants. 

The main conclusions of this work is that to any chemical reaction in a closed system we 
can associate a 2— dimensional equilibrium manifold, and that any chemical reaction can be 
represented as a geodesic in which the initial state is determined by the initial conditions of 
the reaction and the final state corresponds to a coordinate singularity of the thermodynamic 
metric. 
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